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COMPUTATION OF THE LONG RANGE MOTION 
OF A LUNAR SATELLITE 

by 

Isabella J. Cole 
SUMMARY 

Formulae for long periodic perturbations depending on the argument 
of the pericenter are derived by the method of Harmonic Analysis and 
with Elliptic Integrals. 

An approximation of the terms of the Energy Integral that are de- 
pendent on the oblateness was used in the solution with Elliptic Integrals. 
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LIST OF SYMBOLS 


semimajor axis of the satellite's orbit = L 2 //x. 
mean radius of the moon (.1738 decamegametres) . 


eccentricity of the satellite's orbit = 



L /l - e 2 


argument of pericenter of the satellite. 

longitude of the ascending node of the lunar satellite. 

G co s i 

inclination of the satellite's orbit plane to the moon's equatorial plane = 
cos" 1 H/G 

principal part of the oblateness of the moon (2.41 x 10 -4 ). 


3n 2 

c 

8en 



mean anomaly of the lunar satellite. 

'/JL a 

mean motion of the satellite = ^/L 3 . 

mean motion of the moon's mean anomaly = 2.2802713 x io” 3 (radius/centiday) . 
mean motion of the moon's mean longitude = 2.2997150 x 10~ 3 (radians/centidav) . 


vi 



e = ratio of the sum of the masses of the earth and the moon to the mass of the 
earth = 1.0123001. 



(i - gravitational constant times the mass of the moon = 3.6601891 x 10“ 3 (deca- 
megametres 3 / centiday 2 ) . 
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COMPUTATION OF THE LONG RANGE MOTION 
OF A LUNAR SATELLITE 


INTRODUCTION 

In this paper formulae for d-t/dt, dl/dt, -i and h that were derived by: 

1. The method of Harmonic Analysis and 

2. With Elliptic Integrals 
are presented. 


To preserve the continuity of the presentation the formula for g , r) 2 and 
dg/dt from reference 1 are also given. In the solution with Elliptic Integrals 
formulae for d^/dt, dh/dt and dg/dt are presented from reference 3. 


THE HARMONIC ANALYSIS METHOD 

The Harmonic Analysis Method is presented for the computation of the long 
range motion of a Lunar Satellite. The procedure, as outlined below, is well 
adapted to high speed digital computers and attempts to reduce the amount of 
labor involved in the computation. 

a. Derivation of rj 2 as a Cosine Series 


The Hamiltonian, F is given by 


F = F„ 


( 1 ) 


where 


L 3 


F, = n! H 


1 



F 2 = KjL J (5 - 3r/ 2 ) 1 


+ 15(1 -rf) f 1 cos2g .il 1-il UC 


One seeks now an approximation for 1 /rj 3 (dependent on the value of e 2 ) so 
that one has a quadratic in e 2 . The truncated binomial expansion of 1 /if is 
(1 + (3/2)e 2 + (15/8)e 4 ). However, a better approximation for I/17 3 is (1 + (3/2)e 2 
+ 2e 4 ). The following table compares this approximation, the actual value of 
I/77 3 and the binomial expansion of I/77 3 . This approximation was derived 
empirically and is accurate for 0 < e < .5. Figure 1 is a graph of these values. 
Error is defined as actual value/approximate value. 


e 

w 

Author's 

Approximation 

Error 

Binomial 

Expansion 

0 

1 

1 

1 

1 

.1 

1.0151897 

1.0152 

.99998985 

1.0151875 

.2 

1.0631466 

1.0632 

.99994977 

1.063 

.3 

1.1519614 

1.1512 

1.0006614 

1.1501875 

.4 

1.2989160 

1.2912 

1.0059758 

1.288 

.5 

1.5396007 

1.5 

1.0264005 

1.492 


Error 


1.0000022 

1.0001379 

1.0051420 

1.0084752 

1.0319040 


If one approximates I/77 3 by (1 + (3/2)e 2 + 2e 4 ) and substitutes (1 - e 2 ) for 

7? 2 then (neglecting powers of e greater than 4) one has: 

e 4 (3KjL - K 2 (1 - 12 u 2 ) - 15KjL (1 - k 2 ) cos 2g} 

+ e 2 {(KjL + K 2 ) (9v 2 _ 1) + 1 5 Kj L (1 - u 2 ) cos 2g + C} 

+ {2(KjL + K 2 ) (3 v 2 - 1) _C} = 0. (2) 


where C is the value of F 2 with initial conditions. 
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ECCENTRICITY, e 

Figure 1 — Comparison of I/17 3 With Some Approximations of I/77 3 
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C = 6 





3(1 - e 2)3/2 



+ (5/2)K 1 L e 2 0 



(3) 


In equation (2) let 2g take 5 values such that 0 < 2g< n; i.e. for example 0, 
77 / 4 , 77 / 2 , 377/4, 77 . One then has 5 quadratics in e 2 . 

Using the initial conditions of L, G , H and the relationship u = H/L compute 
a , j3 , and y where 

a = coefficient of e 4 term of equation (2) 

fi = coefficient of e 2 term of equation (2) 

y - coefficient of e 2 term of equation (2). 

Solve for rf using the relationship 



and choosing the sign so that 0 < rj 2 < 1. 

Assume that rf can be represented as a cosine series 


i =n 

m. cos 2ig 

i = 0 

i.e. 77 2 = m 0 + nij cos 2g + m 2 cos 4g + m 3 cos 6g 

+ m 4 cos 8 g (for n = 4). (5) 
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Next, construct a 5 x 6 matrix of this shape 


^0 

= m o 

+ 

+ m 2 + m 3 

+ 


i i 2 

s° 

I! 

.5 m - 

.5m 2 _m 3 - 

.5 m. 

4 


i 

~ m 0 


- m 2 

+ m 4 



= m 0 - 

.5 njj - 

.5 m 2 + m 3 - 

-5 m„ 


-nl 

= m 0 


+ m 2 " m 3 

+ m 4 

(6) 


Solve for m 0 . . . m using these relationships. 

I (vl + + y* 3 ) 

"o" -r~ + — r~ 

m i = J< 7 i _7? 4 + ' ^3> 

m 2 = -y 


(t?S - ^ ( 7 ?i 2 - ^ 


m, = 

3 6 


m 4 = t? 2 " m o + m 2 


These values of m 0 . . . m 4 are the coefficients of the cosine series that was 
assumed in equation (5) . 


b. The Angu lar Variables 

Using the Delaunay set of canonical variables (L,G, H, i, g,h) the equations 
of motion for the angular variables become: 
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( 8 ) 


dg _ -dF 
dt " dG 


K, V 


Sv 2 


- 1 | + 5 [1 ) cos 2g 


K„ 


w 


1 - 


Sv 2 


It should be noted that there is the possibility that in certain instances 




Sv 2 


- 1 


1 - — ) cos 2g 

v 4 


K„ 


w 


1 - 


5v 2 


so that dg/dt will be zero and dt/dg will be undefined. 


d l 

dt 


dF 

dL 


K 2 1 




i ±_ 

L 3 


+ 3 K, 


-K^ 2 - 5Kj (2 + v 2 ) cos 2g + 


+ 5K t f cos 2g + - K2 V • J_ + 10K v 2 

L _5 1 2 


- 10 v 2 K, •— 


(9) 


dh 

dt = " dH 


2K 2 H\ i /5K t H\ j_ 

L 2 ) V 5 \ L / T7 2 


r 3K H \ 5K, H ~ 5K, H . 

+ i — i— _ n * + — s 2g — - cos 2 & 

L c + L 7 L 


( 10 ) 


c. g as a function of time, t. 

Using the values of rf obtained from equation (2) compute from equation (8) 

dt J_ 
dg dg 
dt 
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Assume that dt/dg can be represented as a cosine series 


i.e. 


i~n 

p. cos 2ig 

i*0 



dt 

dg 


= P„ + P, cos 2g + P 2 cos 4g + P 3 cos 6g + p 4 cos 8g 


( 11 ) 


Construct a matrix of the same shape and size as the one described as 
equations (6) under the section Derivation of as a Cosine Series, i.e. let p. = 
m. and solve for p 0 . . . p 4 using the relationships given in equation (7) . 

Integrate dt/dg to get 



sin 2g - 



sin 4g - 



sin 6g 


p 4 

sin 8g 

8p o 


or 



sin 2g + — sin 4g + — sin 6g + 


p 4 1 

— sm 8g 

4 


( 12 ) 


Note that g has the holonomic shape, i.e. g can be expressed in terms of time 
and f(g) asg=g + Sf(g) where g is a function of time, t and § is a small 
constant = -l/2p 0 - Equation (12), therefore, can be solved so as to express g 
in terms of 8 and g . (It is assumed that g , f (g) and their derivatives are con- 
tinuous functions). Rewrite g as 


g = g + q sin 2g + r sin 4g + c sin 6g + d sin 8g. 


(13) 


with 


i = 


p 0 


t 


7 


q = 





Then 


g = g + 




sin 2g 


+ 


r + 


B P 


2 + 6 


sin 4g 


+ 





sin 6g 


+ 



D N \ 
+ 2 + 6 ) 


sin 8g 


(14) 


where 


A = - 2(qr + rc + cd) 


B = 2(q 2 - 2qc - 2rd) 


J = 6q(r - d) 


D = 4(2qc + r 2 ) 
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E = 6 r 2 q - c 2 q - d 2 q - - red + 3JL + qrd| 

P = - 24 -j^qrc+qdc+q 2 r + rc 2 +rd 2 - + I_ + j 
M^-54 | r 2 c + red + 3£_ + d 2 c + qdr + q 2 c + ^ 

N = - 48 (- rq 2 + 2qrc + 2r 2 d + re 2 + 2q 2 d + d 3 + 2dc 2 ) 


d. dl/dt and dh/dt as Cosine Series 

Let g, (equation ( 13 )) be such that 0 < g < n i.e. for example 0 , n/4, n/2 , 
377/4, 77 successively and compute g and cos 2 g using equation ( 14 ). Assume 
that co s 2g can be written as the sum of a cosine series 


i a n 

4> i cos 2ig 

i = 0 


to get 


cos 2g = <P 0 + 4> v cos 2g + 4> 2 cos 4g + <t> 3 cos 6g + <& 4 cos 8g ( 15 ) 

Solve for <p 0 . . . 4 > using the procedure described under Derivation of if as 
a cosine series. Compute I/77 3 , r/ 2 cos 2g, 1 /77 s , I/77 2 , and cos (2g/7? 2 ) using the re- 
sults from equations (2) and ( 14 ). Assume that I/77 3 , 1 /77 s , rj 2 cos 2 g, and 
cos( 2 g/) 7 2 ) can b e written as a cosine series. 

For example: 


_L=k„ + k. + k„ + k. ik 

9 V 1 '4 J <4 

% 

_l_ = k 0 + -5 kj - 5 k 2 -k 3 - .5k 4 


9 


from which 


-L-=k 0 -.5k 1 -.5k !+ k 

% 

— — = k Q - kj + k 2 - k 

% 



k 


4 



k 0 +k 2 


cos 2g 


i”n 

L 

i*0 


t j cos 2i g 



i -n 

L 


k. cos 2i g 




— i 


c. cos 2i g 


1 




1 


p. cos 2 ig 


i-n 



rf cos 2 g = 


i-n 

L 


Q { cos 2 ig 


Now if one takes 


d^e 

dt 


<t> • — + r + n T? 2 + A cos 2g + H ' 


„ 1 


and 


with 


+ ¥ C ° S — - 511 rj 2 cos 2g - V • £ T 
V 2 v 2 J 


dh 



•-L- cos 2 g 

V 2 


+ A + 0 cos 2 g 



n = - Kj 

A = - 5Kj (2 + v 2 ) 

3K 2 p 2 
5 = L 


(16) 


( 17 ) 
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W = 10 K t V 2 


r = - 


0 = - 


2K 2 H 

L 2 


5 KjH 


then 


A = 


3 K, H \ 


d i 
d t 


$ ^ ^ o-. cos 2ig + Y 

i e 0 


i^n i= n 

+ / m. cos 2ig + A 


i=0 


4>i cos 2ig 

i^O 


i-n i=n 

+ 3 > P* cos 2ig + W ^ *. cos 2ig 

i = 0 


i -n 

’ E •* 

i *0 


i g n 

- 5n E *■ 

i=0 


cos 2ig 


1 -n 

-* L k ‘ 

i = 0 


cos 2i g 


(18) 


^ r 

dt 


i 3 n 


y-'* P L cos 2ig +0 k. cos 2ig 

i=0 i^O 


+ A + 0 


~ 11 i-n 

y ’ 4> i cos 2i g - © y ' <. cos 2ig 

i=o iTo 


(19) 


Substituting the series in equations (18) and (19) yields the following expres- 
sions for d^/dt anddh/dt 
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w 0 = (Oo- 0 + n (m 0 - 50 o ) + A<£ 0 +Hp 0 + 9(i 0 -k 0 )) + T 


W t = (<&cr. tnCnij - 50 o ) + A<£ 0 + H p, + V( tj -k^) 

W 2 = (<&<r 2 +n(m 2 - 5 0 2 ) + A ^ 2 + Sp 2 + V(i 2 - k 2 )) 

w 3 = (Ocr 3 +n(m 3 - 5 e 3 ) + A0 3 + gp 3 + V( l3 - k 3 )) 

w 4 = (^cr 4 +n (m 4 - 50 4 ) + A0 4 + aP 4 + ^(‘ 4 - k 4 )) 

V 0 =(r^ o+ ©(k 0+ ^ 0 _ t0 ) + A 

v 2 = (r^ + ©ckj + 4> x - < x )) 

v 2 = (V p 2 + ©(k 2 +^ 2 - t 2 )) 

V 3 = (Pp 3 +0(k 3 + <j> 3 - « 3 )) 
v 4 = <T> 4 + 0 (k 4 + <£ 4 - t 4 » 

e. Derivation of i and h 

l and h are obtained by integrating equations (20) and (21) using these 
relationships 

J f(5)d, =flf dS 

dt 


let 


dg _ J_ 

dt Pa 


then 


or 


and 


-Jffl 


dg 


dg 


* - W oo t + p, 


W, 


K - - W 2 

b | 2 " Sln 2e + -Q- slr > 4g 

.) 


■3 _ W 4 

+ — sin 6g + — sin 8g 


h = v oo fc 


W 


00 " W oP 0 


K 

+ p 0 1y sin 2g + 

V 2 . _ V 3 V 8 X 

+ — sm 4g + — sin 6g + — sin 8gJ 


(22) 


(23) 


r oo 


V0P0 


THE SOLUTION WITH ELLIPTIC INTEGRALS 


From equation (1) and the definition of the Delaunay variables one has 


dr? _ 1 dG _ 1 3F „ ,, 2\ 

dt “ L dj - L Si = - 5K l (1 - * ' 

and from equation (3) with H 2 /G 2 = v 2 /-q 2 


v 2 \ . 

sin 

W 


2g 


(24) 


14 





- V 


+ 15 K 1 LC 1 -t, 2 ) (l -iLj cos 2 g ] 


(25) 


from which one finds 


cos 2 g = 


U. 

J 3 KjL V 


¥)*->•■> (-£)•* 


15(1-7,2) 1 


7,2 / 


Substituting this relationship in equation (24) and using 


sin 2 g = + v^l - cos 2 2 g 


(26) 


d V 2 _ dT , 2 dT, _ dT, 


dt d 7 ? dt 


dt 


one has 


dT, 2 

dt 


= 2 7 , 7 , = ± 10 Kj (1 - 7 , 2 ) ( t , 2 _ V 2 ) - /l - cos 2 2 g 


T? 


( 27 ) 


now letting t, 2 = (values of t, 2 for g = 0) = 7 , 2 


K. L | 2 K„ 

I j2 _ 2 > 


c = <(5 - 3 7]*) (3v 2 - r? 2 ) + (3i/ 2 - 7]J) 




K , L ^ 


+ 15 (1 - 7 , 2 ) ( 7,2 - V 2 ) 


( 28 ) 


and 


-jr = + -T ~ j 225(1 - 7 t 2 ) 2 (t , 2 - v 2 ) 2 - 
dt 3 t, y 


2 K 


l K , L 


7 + (5-3TJ 2 ) (tj 2 -3v 2 ) +£2 


-1 2 


vJ 

( 29 ) 


15 



when 1 1 if is approximated by (1 + (3/2)e 2 + 2e 4 ) and equation (29) is expanded 
(neglecting powers of e greater than 8) using the relationships 


+ _ 
dt 3 


77 2 = 

1 - e 2 


d?] 2 

de 2 

de 2 

de 2 

d t 

dt" 

V 

+ A 3 e 6 

+ A 2 e 4 + Aje 2 + A 0 

(1 - e 2 ) 


4 = 225 -^i(l-12, 2 )-3j + - |(1 -9v 2 ) ^_ + l 


2 8K„ 


,, ( K 2 


A = <450 (v 2 - 1) - 2 -4(1 - 3v 2 ) 


2C \ / 2K 2^ 

KjL ) yKjLy 


♦ (rr* 1 ) 


A 2 = <225 (1 - V Z Y - 


4(1 -3v 2 ) 


\ ♦») ^ i i 




A, =-2 2(l-3<) ^ + .) ♦ £ 


(1 -9v 2 ) -4 +1 


A„ = - 2 2(1 - 3y 2 ) 




_c_ 



Let e 2 = y and using the binomial expansion of 

— - — = 1 + y + y 2 + y 3 + y 4 
1 — e 2 

dy/dt becomes (neglecting powers of y greater than 4) 


2k 

^7 = ± ~T /A ;y 4 + A > 3 + Ajy 2 + A Jy + A o 


dt 


( 31 ) 


with 


A 4 - ( A 4 + ^3 + + Aj + A 0 ) 

A 3 = < A 3 + A 2 + A l + V 

A' = (A 2 + Aj + A 0 ) 

A i = (Aj + A 0 ) 


A o = A o 


Equation (31) can be written as: 


•*.± 2 -hvv t 

dt 3 4 


A' A' A' A' 

4 j 3 "2 2 1 A 0 

y + — y J + — y^ + — y+ — 

A A' A' A' 

A 4 A 4 A 4 A 4 


(32) 


and factored as 


2K. 


^7 = ± ~ '/(yj - y) (y 2 - y) (y 3 _ y) (y 4 - y) 


(33) 


where y x » y 2 » y 3 > y 4 are the zeros of the quartic in equation (32) from which 


. 2K - /A i dt 

^ 3 
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Q r (yj - y) (y 2 - y) (y 3 - y) (y 4 - y) 

such that 

y 4 > y 3 > y 2 > Vi > y 0 - 

let 

2 Kj /a; 

* 3 t 

k2 = (y 4 ~ Vi) (y 3 ~ y 2 ) 

(y 4 - y 2 ) (y 3 - y^ 

, (y 4 - y 2 ) (yi - y) 

sn z u = — 

(y 4 - y x ) (y 2 - y) 

2 

g - — - — - 

^(y 4 - y 2 ) (y 3 - yj) 

and 


I 


Vl dy 

/O' 


gF(0, k) 


however one has 



now let 


so 


and 



dy 

/Q _ 


(u ± V) 



dy 


(u ± v) = sn' 1 (y 2 ) - sn' 1 (y) 
sn (u ± v) = sn (sn' 1 y 2 ) - sn (sn -1 (y)) 

sn (u ± v) = y 2 - y or y = y 2 - sn (u + v) 


y = y 2 - 


(sn ucnvdnv ± snvcnudnu) 
1 - k 2 sn 2 u sn 2 v 


(35) 


now 


y = e 2 = 1 - T} 2 and r} 2 - 1 - y 


or 


v 2 = l - y 2 + 


sn u cn vdn v ± sn v cn u dn u 
1 - k 2 sn 2 u sn 2 v 


(36) 


Once rj 1 is known, dg/dt, dh/dt and d'C/dt can be integrated us ing rf as the 
independent variable. For purposes of continuity with reference 3 one calls rj 2 = 
x then 



j* * dx 

Jx 0 (1 - X) v/x Q(x) 


19 


where 


4G ' f 

j X 


dx 


+ G, 


+ G, 


1 


G <1 


4 


(X 

- V 2 ) /x Q(x) 


dx 

/x 

Q(x) 


dx 

(1 

- x) /Q(x) 


dx 

(X 

- v 2 ) /Q(x) 


dx 


+ G, 


+ G, 


x„ x /0(F) 
dx 


j 


x 0 F /0(X) 


f 


dx 


/Q(x) 


12 Kj /6* 


G o = [j + 2 (1-3^)^ 

Gj = -^+2v J (3^-5)^] 

G 2 = - [y + 2 C 1 + 3v 2 ) Kj 


2 K~ 


G 3 = — (l-3v 2 ) 


(37) 
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with 


4K> 

G 4 = — - 

4 L v 2 

G = — — ( al ~ 3v 4 - 2^ 

5 " L V - 2 > 

-K, 

G 6 = — (7 + 6v 2 ) 

9v 2 K, 
g 7 =-ir- 


M = L 0 t + Jl, 



dx 


0 (1 - X) /QOO 


+ L- 


I 


dx 


* 0 x /* Q(x) 


+ L, 


I' 


dx 


v^Tqoo 


+ u 




dx 


(1 - x) ^TqooJ \ 4K, /6~ 


L i 


L- = 


( 2 „ C \ 

In + — K, 

\ 3 1 3L / 

Kj (2 - 6v 2 ) c 




3L 


K, 

--(7-6^) 


L, - + 


7 v 2 K„ 


L a = - 


2 K 2 (1 - 3v 2 ) 
3L 
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0 


Ah = H 0 t 



with 


where 



Q(x) = (x - x,) (x - x 2 ) (x - x 3 ) (x - x 4 ) 

The integrals needed are of two types: 

I. Integrals involving / Q(x) 


for 

and 


1 . 


2 . 


f 

/ 


dx 

(1 - X) / Q(x) 


dx 

(^ 2 - X) /ooo 


see reference 5 page 102 # 251.39 



— - dt 

'y (P - ty /(a-t) (b-t) (c-t) (t-d) 

g f Ul (1-a 2 sn 2 u) m du _ 
(p-d) m J o (l-a 2 sn 2 u) m 


with m = 1 


see #340.01 page 205 reference 5 


with 


g 

a 3 (P" d ) 



a 2 ) II (<p, a 2 , k) + a 2 uj 


sn 2 u = ( a - c > ( d -0 
(a -d) (c-t) 


g= 


a 2 = 


/(a - c) (b - d) 
(p - c) (a - d) 


(p - d) (a - c) 


; p t d, a j- c 


cp - am 


Uj = sin" 1 


(a - c) (d - y) 
(a - d) (c - y) 


2 _ a — d 


a* - 


a - c 


, u = F (q>, k) 


t = d 


1 - a 2 sn 2 u\ 

V 1 - a 2 sn 2 u/ 
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for 


3. 


4. 


5. 


I 

I 

S 


dx 



dx 

X 2 / Q(x) 


dx 

X 3 /Q(x) 


see reference 5 page 99 No. 251.04 


f 


dt 


t m /(a - t) (h - t) (c - t) (d - t) 


with 


g 

d m 


f 

•'o 


(1 - a 2 sn 2 u) m 
(1 - a 2 sn 2 u) n 


du 


m = 1, 2, 3. 


2 ca' 

ai - 

2 A 


with m = 1 


t = d 


^1 = a 2 sn 2 u' 
l 1 - a 2 sn 2 u i 


see #340.01 page 205 reference 5 


a 2 d 


[( a 2 " °- 2 ) n (< P . a 2 ’ k ) + a2 ^ 


24 



with m = 2 


see #340.02 page 205 reference 5 

~ [* 4 u + 2a2 (a 2 - a 2 ) V, + (a 2 - a*) 2 vj 

where 

V x = ri(q>, a 2 , k) 

#336.01 page 201 reference 5 


V 


2 


1 

2 ( a 2 - !) ( k2 - a 2 ) 

+ (2 a 2 k 2 + 2 a 2 - 


|a 2 E(u) + (k 2 - a 2 ) u 
a\ - 3 k 2 ) II (<p, a 2 k) 


a* sn u cn u dn u 
1 - a 2 sn 2 u 

#336.02 page 201 reference 5 
with m = 3 

see #340.03 page 205 reference 5 


g 

af d 3 


a 6 u + 3 a 4 (a 2 - a 2 ) V x + 3 a 2 


( a 2 


-a 2 ) 


2 ' 2 V„ 


+ ( a 



where 


V 


3 


g 

4 (1 - a 2 ) (k 2 - a 2 ) d 3 



2 (a 2 k 2 + a 2 - 3 k 2 ) Vj 
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+ 3 (a* - 2 a 2 k 2 - 2 a 2 + 3 k 2 ) V 2 + 




For these integrals one fits, by the method of least squares, a polynomial 
of degree 2 to the expression for x Q(x), (0 < x < x.). The procedure is 



z = (y + *!> (y + a 2 ) (y + a 3 ) (y + a 4 ) (y + a 5 ) a i = i 



by using the transformation x = y + 1 with - 1 < y < - a 5 

2. Assign y values such that 

y i = (- 1 + i o i i < io. 

3. Compute the normal matrix 


11 a' 

M 

cr 

2 y 2 c' 

IZ 

2y a' 

2 y 2 b' 

— y 3 c' 

ly Z 

I y 2 a' 

Sy 3 b' 

£ y 4 c' 

2y 2 Z 


4. Solve this system for a', b' c'. These values are the coefficients of the 
approximating least squares polynomial 

t it *2 

a + b y + c y z 
or 

a + b x + c x 2 

with 

a = a' - b' + c' 
b = b' - 2 c' 
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Set a + bx + cx 2 = X and one has that 


6. = 

r dx 

dx 


' (1 - x) / x Q(x) 

; (i - x> /x 

7. = 1 

dx 

r dx 

j 

(x - V 2 ./x Q(x) 

^ (X - Z^ 2 ) /X 

8. = | 

dx ( 

♦ 

cx 

J 

/xQ(x) J 

X 

9. = 

dx / 

dx 


x/xQ(x) J 

X /X 

10. = 1 

dx __ j" 

dx 

J 

x 2 A Q(x) J 

x 2 /X 


for 6 and 7. 

Use #195, 196 of reference 6 page 27 

ifk*0 J log (2k +/3v- 2 b' /kX) 

J v /X /k v 


or 



tan 1 


2 k + ft v 
2b' /-kX 


sin 

v^k 


2k + yg v 

b' v /-q 
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V = (a* + b' x) 


/S = bb' - 2a' c 
k = ab' 2 - a' bb' + c a' 2 
q = 4 ac - b 2 


if k = 0 


dx - 2 b' ✓' X 


J-V 


for 8 see #160, 161 of reference 6 page 23 


if c > 0 


I 


-fo = _L log f/X + x/c + 

V X V 2 /c. 


if c = 0 


/■ 


dx _ 2 . X 

7x ” “ 


if c < 0 


r dx 

J 7f 


dx “1 . _ i [ 2 c x + b\ 

- sm 1 [ 1 


■/ -c 


v ~ a 


Isinh-' 2CX+ b 


VC 


✓ q 


for 9. See #182, 183 of reference 6 pages 25, 26. 


if a > o 


[* dx 1 / /X + ✓ a b 

j 777” " 77 ° g V * + 777 
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if a = 0 


if a < 0 


J 

J 


dx 

x /X 


dx 
x /X 



1 

/^i 


sin 1 


^bx + 2 a\ 

wry 


or 


L sinh' 1 

/I 


/2 a + bx\ 

\~ 7q j 


for 10. See #186 of reference 6 page 26. 


dx ~ /X b f dx 

J x 2 7x ax 2 a J x/l' 

An example of the fitting procedure is as follows: 


let 


such that 


Xj = .53 x 3 = .23 

x 2 = .31 x 4 = .12 

a = 1 c = .69 

b = .47 d = .77 

e = .88 


y = -1, -.988, -.976 • • • -.88 

and the normal matrix is 


11a' 

-10.34b' 

9.73544c' 

4.8922239 x 

10' 4 

-10.34a' 

9.73544b' 

-9.1810928c' 

-4.6383256 x 

*r 

i 

O 

i— i 

9.73544a' 

-9. 1810928b' 

8.6722565c' 

4.4007739 x 

10' 4 
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whose solution gives as the approximate formula 


Z - (178.02728 -378.13608 y -199.80525 y 2 ) 10' 4 
or 

Z ~ (.30355 + 21 . 474220 x -199.80525 z 2 ) 10' 4 

The original value of Z and those obtained by substitution in the approximating 
formula are plotted below. 



- 1.000 - 0.976 - 0 952 - 0.928 - 0.904 - 0.880 


V 


31 




REFERENCES AND BIBLIOGRAPHY 


1. Frost, Frances A., Long Term Motion of a Lunar Satellite, GSFC X-547-65-388 
October 1965. 

2. Brown, E. W., "The Stellar Problem of Three Bodies", Monthly notices of the 
Royal Astronomical Society Vol. 97, 56-66, 116-127, 388-395 (1936). 

3. Giacaglia, Georgio E. O., Murphy, James P., and Felsentreger, Theodore L. 
in Colloboration with Standish, E., Myles, Jr. and Velez, Carmelo E., The 
Motion of a Satellite of the Moon , GSFC X-547-65-218, June 1965. 

4. Kozai, Y., "Motion of a Lunar Orbiter", Journal of the Astronomical Society 
of Japan, Vol. 15, 301-31 (1963). 

5. Byrd, P., and Friedman, M., Handboo k of Elliptic Integrals for Engineers and 
Physicists , Springer- Verlag, 1954. 

6. A Short Table of Integrals, by B. O. Peirce, Ginn and Co. 1910. 


32 





